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Solar neutrinos from the decay of *B have been detected at the Sudbury Neutrino Observatory 
(SNO) via the charged current (CC) reaction on deuterium and by the elastic scattering (ES) of elec- 
trons. The CC reaction is sensitive exclusively to i^e's, while the ES reaction also has a small sensi- 
tivity to i/^'s and Vt^s. The flux of i^e's from *B decay measured by the CC reaction rate is (fP'^{ve) — 
1.75 ±0.07 (stat.)t°;i^ (sys.)±0.05 (theor.) x 10"^ cm s . Assuming no flavor transformation, the 
flux inferred from the ES reaction rate is (p^^if^) = 2.39 ± 0.34 (stat.)+° (sys.) x lO" cm"^s"^ 
Comparison of (fP'^{ve) to the Super-Kamiokande Collaboration's precision value oi<f)^^{vx) yields a 
3.3(7 difference, assuming the systematic uncertainties are normally distributed, providing evidence 
that there is a non-electron flavor active neutrino component in the solar flux. The total flux of 
active *B neutrinos is thus determined to be 5.44 ±0.99 x 10® cm~'^s~^, in close agreement with the 
predictions of solar models. 
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Solar neutrino experiments over the past 30 years 
have measured fewer neutrinos than are predicted by 
models of the Sun . One explanation for the deficit is 
the transformation of the Sun's electron-type neutrinos 
into other active flavors. The Sudbury Neutrino Obser- 
vatory (SNO) measures the ®B solar neutrinos through 
the reactions: 



TABLE I. Data reduction steps. 



Analysis step 



Number of events 



i^e + d ' 
Vx + d 
Vx + 



p -l-p -I- e 
Vx + e~ 



(CC) 
(NC) 
(ES) 



The charged current reaction (CC) is sensitive exclu- 
sively to electron-type neutrinos, while the neutral cur- 
rent (NC) is sensitive to all active neutrino flavors {x = 
e,fi,T). The elastic scattering (ES) reaction is sensitive 
to all flavors as well, but with reduced sensitivity to z/^ 
and Vt- By itself, the ES reaction cannot provide a mea- 
sure of the total flux or its flavor content. Comparison 
of the *B flux deduced from the ES reaction assuming no 
neutrino oscillations {(p^^{i'x)), to that measured by the 
CC reaction (0*^*-^(i'e)) can provide clear evidence of fla- 
vor transformation without reference to solar model flux 
calculations. If neutrinos from the Sun change into other 
active flavors, then 0^^(z^e) < 0^^(i/a;). 

This Letter presents the first results from SNO on the 
ES and CC reactions. SNO's measurement of (f>^^{i'x) 
is consistent with previous measurements described in 
Ref Q. The measurement of cf)'^'^ (ve) , however, is sig- 
nificantly smaller and is therefore inconsistent with the 
null hypothesis that all observed solar neutrinos are i^e- 
A measurement using the NC reaction, which has equal 
sensitivity to all neutrino flavors, will be reported in a 
future publication. 

SNO is an imaging water Cerenkov detector lo- 
cated at a depth of 6010 m of water equivalent in the 
INCO, Ltd. Creighton mine near Sudbury, Ontario. It 
features 1000 metric tons of ultra-pure D2O contained 
in a 12 m diameter spherical acrylic vessel. This sphere 
is surrounded by a shield of ultra-pure II2O contained 
in a 34 m high barrel-shaped cavity of maximum diame- 
ter 22 m. A stainless steel structure 17.8 m in diameter 
supports 9456 20-cm photomultiplier tubes (PMTs) with 
light concentrators. Approximately 55% of the light pro- 
duced within 7 m of the center of the detector will strike 
a PMT. 

The data reported here were recorded between Nov. 2, 
1999 and Jan. 15, 2001 and correspond to a live time of 
240.95 days. Events are defined by a multiplicity trig- 
ger of 18 or more PMTs exceeding a threshold of ^ 0.25 
photo-electrons within a time window of 93 ns. The trig- 
ger reaches 100% efficiency at 23 PMTs. The total in- 
stantaneous trigger rate is 15-18 Hz, of which 6-8 Hz is 
the data trigger. For every event trigger, the time and 
charge responses of each participating PMT are recorded. 

The data were partitioned into two sets, with approx- 
imately 70% used to establish the data analysis proce- 



Total event triggers 
Neutrino data triggers 
A'hit >30 

Instrumental background cuts 
Muon followers 
High level cuts'* 
Fiducial volume cut 
Threshold cut 



355 320 964 
143 756 178 
6 372 899 
1 842 491 
1 809 979 
923 717 
17 884 
1 169 



Total events 



1 169 



''Reconstruction figures of merit, prompt light, and {9ij}- 

dures and 30% reserved for a blind test of statistical bias 
in the analysis. The analysis procedures were frozen be- 
fore the blind data set was analyzed, and no statistically 
significant differences in the data sets were found. We 
present here the analysis of the combined data sets. 

Calibration of the PMT time and charge pedestals, 
slopes, offsets, charge vs. time dependencies, and second 
order rate dependencies are performed using electronic 
pulsers and pulsed light sources. Optical calibration is 
obtained using a diffuse source of pulsed laser light at 
337, 365, 386, 420, 500 and 620 nm. The absolute energy 
scale and uncertainties are established with a triggered 
^^N source (predominantly 6.13-MeV 7's) deployed over 
two planar grids within the D2O and a linear grid in 
the H2O. The resulting Monte Carlo predictions of de- 
tector response are tested using a ^^^Cf neutron source, 
which provides an extended distribution of 6.25-MeV 7 
rays from neutron capture, and a ■^H(p, 7)^He |l^] source 
providing 19.8-MeV 7 rays. The volume- weighted mean 
response is approximately nine PMT hits per MeV of 
electron energy. 

Table || details the steps in data reduction. The first 
of these is the elimination of instrumental backgrounds. 
Electrical pickup may produce false PMT hits, while elec- 
trical discharges in the PMTs or insulating detector ma- 
terials produce light. These backgrounds have character- 
istics very different from Cerenkov light, and are elimi- 
nated using cuts based only on the PMT positions, the 
PMT time and charge data, event-to-event time correla- 
tions, and veto PMTs. This step in the data reduction 
is verified by comparing results from two independent 
background rejection analyses. 

For events passing the first stage, the calibrated times 
and positions of the hit PMTs are used to reconstruct the 
vertex position and the direction of the particle. The re- 
construction accuracy and resolution are measured using 
Compton electrons from the ^^N source, and the energy 
and source variation of reconstruction are checked with 
a *Li f3 source. Angular resolution is measured using 
Compton electrons produced more than 150 cm from the 
^^N source. At these energies, the vertex resolution is 
16 cm and the angular resolution is 26.7 degrees. 
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An effective kinetic energy, Tes, is assigned to eacfi 
event passing the reconstruction stage. T^g is calculated 
using prompt (unscattered) Cerenkov photons and the 
position and direction of the event. The derived energy 
response of the detector can be characterized by a Gaus- 
sian: 



R{EcS, Ee) — 



1 



2T:aE{Ee) 



exp[ 



'2^ asiE,) 



■f 



where Ee is the total electron energy, EcS — Toff + me, 
imdaE^E^) = (-0.4620 + 0.5470\/eI+0.008722£;e) MeV 
is the energy resolution. The uncertainty on the energy 
scale is found to be ±1.4%, which results in a flux un- 
certainty nearly 4 times larger. For validation, a second 
energy estimator counts all PMTs hit in each event, iVhit, 
without position and direction corrections. 

Further instrumental background rejection is obtained 
using reconstruction figures of merit, PMT time resid- 
uals, and the average angle between hit PMTs ((6*^)), 
measured from the reconstructed vertex. These cuts test 
the hypothesis that each event has the characteristics 
of single electron Cerenkov light. The effects of these 
and the rest of the instrumental background removal cuts 
on neutrino signals are quantified using the *Li and ^^N 
sources deployed throughout the detector. The volume- 
weighted neutrino signal loss is measured to be 1.4lo;g% 
and the residual instrumental contamination for the data 
set within the D2O is < 0.2%. Lastly, cosmic ray induced 
neutrons and spallation products are removed using a 20 s 
coincidence window with the parent muon. 

Figure |l| shows the radial distribution of all remaining 
events above a threshold of rcff>6.75 MeV. The distri- 
bution is expressed as a function of the volume- weighted 
radial variable {R/ Rav)^ , where Raw = 6.00 m is the 
radius of the acrylic vessel. Above this energy threshold, 
there are contributions from CC events in the D2O, ES 
events in the D2O and H2O, a residual tail of neutron 
capture events, and high energy 7 rays from radioactiv- 
ity in the outer detector. The data show a clear signal 
within the D2O volume. For (R/Rav)^ > 1-0 the distri- 
bution rises into the H2O region until it is cut off by the 
acceptance of the PMT light collectors at i? 7.0 m. A 
fiducial volume cut is applied at i? = 5.50 m to reduce 
backgrounds from regions exterior to the D2O, and to 
minimize systematic uncertainties associated with optics 
and reconstruction near the acrylic vessel. 

Possible backgrounds from radioactivity in the D2O 
and H2O are measured by regular low level radio- assays 
of U and Th decay chain products in these regions. The 
Cerenkov light character of D2O and H2O radioactiv- 
ity backgrounds is used in situ to monitor backgrounds 
between radio-assays. Low energy radioactivity back- 
grounds are removed by the high threshold imposed, as 
are most neutron capture events. Monte Carlo calcula- 
tions predict that the H2O shield effectively reduces con- 
tributions of low energy (< 4 MeV) 7 rays from the PMT 




FIG. 1. Distribution of event candidates with rcff>6.75 
MeV as a function of the volume weighted radial variable 
(R/Rav)"^. The Monte Carlo simulation of the signals, 
weighted by the results from the signal extraction, is shown 
as a histogram. The dotted line indicates the fiducial volume 
cut used in this analysis. 



array, and these predictions are verified by deploying an 
encapsulated Th source in the vicinity of the PMT sup- 
port sphere. High energy 7 rays from the cavity are also 
attenuated by the H2O shield. A limit on their leakage 
into the fiducial volume is estimated by deploying the 
^^N source near the edge of the detector's active volume. 
The total contribution from all radioactivity in the de- 
tector is found to be <0.2% for low energy backgrounds 
and <0.8% for high energy backgrounds. 

The final data set contains 1169 events after the fidu- 
cial volume and kinetic energy threshold cuts. Figure |^ 
(a) displays the distribution of cos 6q , the angle between 
the reconstructed direction of the event and the instanta- 
neous direction from the Sun to the Earth. The forward 
peak in this distribution arises from the kinematics of 
the ES reaction, while CC electrons are expected to have 
a distribution which is (1 — 0.340 cos 6*0 ) yj], before ac- 
counting for detector response. 

The data are resolved into contributions from CC, ES, 
and neutron events above threshold using probability 
density functions (pdfs) in roff7 cos0q, and (i?/i?Av)'^, 
generated from Monte Carlo simulations assuming no fla- 
vor transformation and the shape of the standard 
spectrum [hep neutrinos are not included in the 
fit). The extended maximum likelihood method used 
in the signal extraction yields 975.4±39.7 CC events, 
106.1±15.2 ES events, and 87.5±24.7 neutron events for 
the fiducial volume and the threshold chosen, where the 
uncertainties given are statistical only. The dominant 
sources of systematic uncertainty in this signal extrac- 
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TABLE II. Systematic error on fluxes. 
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FIG. 2. Distributions of (a) cos 6^0, and (b) extracted ki- 
netic energy spectrum for CC events with R <5.50 m and 
rcff>6.75 MeV. The Monte Carlo simulations for an undis- 
torted *B spectrum are shown as histograms. The ratio of 
the data to the expected kinetic energy distribution with cor- 
related systematic errors is shown in (c). The uncertainties 
in the *B spectrum jl^ have not been included. 



tion are the energy scale uncertainty and reconstruction 
accuracy, as shown in Table The CC and ES signal de- 
composition gives consistent results when used with the 
iVhit energy estimator, as well as with different choices 
of the analysis threshold and the fiducial volume up to 
6.20 m with backgrounds characterized by pdfs. 

The CC spectrum can be extracted from the data by 
removing the constraint on the shape of the CC pdf and 
repeating the signal extraction. 

Figure |^ (b) shows the kinetic energy spectrum with 
statistical error bars, with the spectrum of Ortiz et 
al. scaled to the data. The ratio of the data to the 
prediction |^] is shown in Figure || (c). The bands rep- 
resent the la uncertainties derived from the most signif- 
icant energy-dependent systematic errors. There is no 
evidence for a deviation of the spectral shape from the 
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predicted shape under the non-oscillation hypothesis. 

Normalized to the integrated rates above the kinetic 
energy threshold of Tes = 6.75 MeV, the measured 
®B neutrino fluxes assuming the standard spectrum 
shape are: 



1.75 ±0.07 (stat. 
xlO*^ 



tail (sys.) ± 0.05 (theor.) 



cm 



2.39±0.34(stat.)t[3:i^ (sys.) x 10^ cm^^g 



where the theoretical uncertainty is the CC cross section 
uncertainty [l^ . Radiative corrections have not been 
applied to the CC cross section, but they are expected 
to decrease the measured {i/e) flux Q by up to a 
few percent. The difference between the flux deduced 
from the ES rate and that deduced from the CC rate in 
SNO is 0.64 ± 0.40 x 10^ cm'^s'^, or 1.6ct. SNO's ES 
rate measurement is consistent with the precision mea- 
surement by the Super-Kamiokande Collaboration of the 
^B flux using the same ES reaction ||]: 



/,fS(j,^) = 2.32 ±0.03 (stat.)t°:g? (sys.) 



X 10^ 



The difference between the flux (p^^li'x) measured by 
Super-Kamiokande via the ES reaction and the (f)^'-^ (vg) 
flux measured by SNO via the CC reaction is 0.57 ± 
0.17 X 10^ cm^^s"^, or 3.3ct ||1^, assuming that the sys- 
tematic errors are normally distributed. The probability 
that a downward fluctuation of the Super-Kamiokande 
result would produce a SNO result > 3.3ct is 0.04%. For 
reference, the ratio of the SNO CC ^B flux to that of the 
BPBOl solar model § is 0.347±0.029, where all uncer- 
tainties are added in quadrature. 

If oscillation solely to a sterile neutrino is occur- 
ring, the SNO CC-derived ^B flux above a threshold of 
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6.75 MeV will be consistent with the integrated Super- 
Kamiokande ES-derived *B flux above a threshold of 8.5 
MeV Adjusting the ES threshold g this derived 

flux difference is 0.53±0.17 x 10^ cm-'^s-^, or 3. la. The 
probability of a downward fluctuation > 3.1a is 0.13%. 
These data are therefore evidence of a non-electron active 
flavor component in the solar neutrino flux. These data 
are also inconsistent with the "Just-So^" parameters for 
neutrino oscillation | p^ . 

Figure ^ displays the inferred flux of non-electron fla- 
vor active neutrinos {4>{i'^r)) against the flux of electron 
neutrinos. The two data bands represent the one stan- 
dard deviation measurements of the SNO CC rate and 
the Super-Kamiokande ES rate. The error ellipses rep- 
resent the 68%, 95%, and 99% joint probability contours 
for (jyivp) and (tiiv^r)- The best fit to Mvut) is: 



E 

4- 



(j){v^r) = 3.69 ± 1.13 X l(f cm^^s^^ 
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FIG. 3. Flux of *B solar neutrinos which are ^ or r flavor 
vs. the flux of electron neutrinos as deduced from the SNO 
and Super-Kamiokande data. The diagonal bands show the 
total *B flux <j){vx) as predicted by BPBOl (dashed lines) and 
that derived from the SNO and Super-Kamiokande measure- 
ments (solid lines). The intercepts of these bands with the 
axes represent the ilcr errors. 



The total flux of active °B neutrinos is determined to 



be: 



5.44 ±0.99 X 10" cm^^s 



This result is displayed as a diagonal band in Fig. and 
is in excellent agreement with predictions of standard 
solar models . 

Assuming that the oscillation of massive neutrinos 
explains both the evidence for electron neutrino flavor 
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change presented here and the atmospheric neutrino data 
of the Super-Kamiokande collaboration |l^ , two separate 
splittings of the squares of the neutrino mass eigenvalues 
are indicated: < 10~'^eV^ for the solar sector [ p"7|jT^ ] and 
~ 3.5 X 10~^eV^ for atmospheric neutrinos. These re- 
sults, together with the beta spectrum of tritium [ p2[ , 
limit the sum of mass eigenvalues of active neutrinos to 
be between 0.05 and 8.4 eV, corresponding to a constraint 
of 0.001 < rijy < 0.18 for the contribution to the critical 
density of the Universe |2^,^ . 

In summary, the results presented here are the first 
direct indication of a non-electron flavor component in 
the solar neutrino flux, and enable the flrst determination 
of the total flux of neutrinos generated by the Sun. 
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